New results for the virial coefficients of D— dimensional hard spheres 
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Exact results are given for the fourth virial coefficient of hard spheres in even dimensions up 
through 12. The fifth and sixth virial coefficients are numerically computed for dimensions 2 through 
50 and it is found that the sixth virial coefficient is negative for D > 6. Numerical studies are made 
of the contributing Ree Hoover diagrams up to order 17. It is found for D > 3 that for large order 
a class of diagrams we call "loose packed" dominates and the rate of growth of these diagrams is 
used to study bounds on the radius of convergence. 
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The virial series for the pressure 

oo 

for hard particles of diameter a where the pair poten- 
tial £7(r) is +oo for |r| < a and otherwise, has been 
studied in dimensions two and three for over 100 years. 
However, despite the long history of this problem, the 
only rigorous information of the radius of convergence 
comes from the upper bound of Lebowitz and Penrose [1] 
that \B k /B%~ l \ < 13.8 k - 1 /k which in terms of the 
packing fraction r\ = B2p/2 D ~ 1 gives convergence for 
|ry| < 0.145/2 15 . For D = 3 this gives a lower bound 
of 0.018 which is much less than the packing fraction at 
which freezing occurs of rjf = 0.49 [2, 3]. 

Our numerical knowledge of the virial expansion for 
D = 2, 3 is compactly summarized in Table I. Virial 
coefficients have also been calculated for hard spheres in 
dimensions higher than three in [4-7]. 

TABLE I: Bu/B^ 1 for k = 3, . . . , 8 in D = 2, 3. 

discs spheres 

B 2 tto- 2 /2 2trt 3 / 3 

Bs/B 2 2 f - ^ [8] | [9] 

Bt/B'i 0.5322318 ■ ■ ■ [10, 11] 0.2869495 • • • [9, 12, 13] 

Bs/Bl 0.33355604(4) [14, 15] 0.110252(1) [14, 16] 

B 6 /Bl 0.19883(1) [14, 17] 0.038808(55) [14, 18] 

B 7 /B$ 0.114877(11) [18-20] 0.013046(22) [18-20] 

Bs/Bl 0.065030(31) [18, 20] 0.004164(16) [18, 20] 



There have been many attempts to estimate a radius of 
convergence and to find an approximate low density form 
for the equation of state based on these first few virial co- 
efficients. In three dimensions these approximates may 
be grouped into three classes in terms of the location 
of the leading singularity in the 77 plane. Either there 
is (1) a high order pole [21-24] at 77 = 1, (2) a simple 
pole [3, 14, 18, 25, 26] at or near the packing fraction 



Vcp = 0.74048 • • • of closest packed spheres, or (3) a frac- 
tional power law divergence [27-30] at or near the "ran- 
dom close packed" density rj rcp — 0.64 as defined by [31]. 

All of these approximates have radii of convergence 
greater than the endpoints of the numerically deter- 
mined [2, 3] first order phase transition rjf = 0.49 and 
r)s = 0.54; because of this it is often assumed that the 
virial expansion for hard spheres is analytic at the freez- 
ing density 77/. This analyticity assumption is incorpo- 
rated in most of the phenomenological theories used to 
describe freezing [32-34] . There is clearly a large discrep- 
ancy between the lower bound on the radius of conver- 
gence and the assumption of analyticity at rjf. 

Furthermore all the known virial coefficients for D = 
2, 3, 4, 5 are positive and this positivity is often built into 
the approximate equations of state by having the lead- 
ing singularity on the positive real axis. But the pos- 
sibility of negative virial coefficients was suggested by 
Temperley [35] as far back as 1957, and in 1964 it was 
shown [4] that B4 is negative for D > 8. The sixth and 
seventh virial coefficients are negative for parallel hard 
cubes [36] and oscillatory signs are found in models of 
hard squares [37] and hard hexagons [38, 39]. Conse- 
quently there is no a priori reason to expect the leading 
singularity to be on the real axis. 

We have studied these questions of sign change and 
rate of growth of the virial coefficients by numerically, 
and for B4 analytically, evaluating selected Ree Hoover 
diagrams for orders up through 17. We use the defini- 
tions and notations adapted from [40]. In particular, for 
diagrams of k points which contribute to B k , each point 
is connected to each other point either by an / bond that 
corresponds to the function /(r) which is —1 for |r| < a 
and otherwise, or by a bond / = 1 + /. In drawing 
the diagrams we need only specify either / bonds (repre- 
sented by solid lines) or / bonds, not both. For example 
the fourth virial coefficient is given by 

s 4 = -0-- = -^--n- (2) 

4 8 4^ 8 U V ' 

The first expression in Eq. 2 is the expansion in Ree- 



2 



Hoover graphs with / bonds, while the second is the same 
expansion but with only / bonds shown. 

The integrals in Eq. 2 have been previously evaluated 
analytically only for D = 2 [10, 11] and 3 [9, 12, 13]. 
We have extended these analytic computations for even 
dimensions up through D = 12. The results are shown 
in Table II. Details of the computation will be published 
elsewhere. 

TABLE II: Analytical results for B4/B2 in even dimensions. 



D Analytic Value 



Numerical Value 



2 

4 

6 

8 

10 

12 



2 - 
2 - 
2 - 
2 - 
2 - 
2 - 



832 

Is?? 1 

38848 



9\/3 , 

4tt ' 
8JV3 1 

IOtt ' 
2511V / 3 

280-n- 
2673 y/3 1 

280-zr 1B280657T 5 
2187y / 3 1 11565604768 

2207T t xtrnnum^ 



17605024 
606375-n- 2 
49048616 



0.53223180 
0.15184606 
0.03336314 
-0.00255768 
-0.01096248 
-0.01067028 



To investigate the phenomena of negative virial co- 
efficients further we have made a Monte-Carlo evalua- 
tion of i?4,i?5, and Bq for dimensions up to 50. The 
method used allows the calculation of Bk for dimensions 
D > k — 1, including non- integer dimensions, and will be 
reported elsewhere. The results for integer dimensions 
up to 12 are shown in Table III where the entries for Bq 
in D = 4 and 5 are from [7]. For dimensions higher than 
12 these Bk/B^ 1 approach zero in a monotonic fash- 
ion. From an interpolation of these results we see that 
B4 becomes negative at D = 7.73, B§ becomes negative 
at D ~ 5.3, and that while B^jB\ is always positive it 
is not monotonic. The fact that the zero crossing of Bk 
has decreased from 7.73 to 5.3 as D increases from 4 to 
6 suggests that it is plausible that the zero crossing for 
B 8 occurs close to D = 4. 



TABLE III: Numerical values for B A /Bl, B 5 /B% and B 6 /B%. 
Local minima and maxima are underlined. 



D 


Ba/BI 


B 5 /Bl 


B & /Bl 


3 


0.2869495 ■ • • 


0.110252(1) 


0.03881(6) 


4 


0.1518460- •• 


0.03565(5) 


0.00769(3) 


5 


0.075978(4) 


0.01297(1) 


0.00094(3) 


6 


0.03336314 • • 


• 0.007528(8) 


-0.00176(2) 


7 


0.009873(4) 


0.007071(7) 


-0.00352(2) 


8 


-0.0025576 ■ • • 


0.007429(6) 


-0.00451(2) 


9 


-0.008575(3) 


0.007438(6) 


-0.00478(1) 


10 


-0.0109624 ■ • • 


0.006969(5) 


-0.00452(1) 


11 


-0.011334(3) 


0.006176(4) 


-0.00395(1) 


12 


-0.0109624- •• 


0.005244(4) 


-0.003261(7) 



To proceed further we separate the Ree-Hoover dia- 
grams with k points into classes with to < k points which 
are the end points of / bonds. We designate the value of 



a diagram in this class, including its combinatorial coeffi- 
cient, as Bk[m, i] where i is an arbitrary label specifying 
the graph with given k and to. If the graph B k [to, i] exists 
for k — to then it will continue to exist for k > to. There 
is one graph in the class Bk[0, i], one in the class Bfe[4, i] 
three in the class Bk[5,i], and 18 in the class Bk[6,i\. 

For the diagram with to = the absence of / bonds 
forces the points to all lie within a distance a of each 
other. We define any sequence of diagrams with increas- 
ing number of points as being "close packed" if the max- 
imal volume of the convex hull of k points approaches a 
constant as k — > 00. As such, any sequence with to and 
i fixed is a close packed sequence of diagrams. We have 
numerically studied Bfe[0, 1], £?fc[4, 1], and all graphs in 
Bk[5,i] for many values of k and D. Two examples are 
shown in Tables IV and V. 



TABLE IV: B k [0. 



0/B2 



D 



D ■■ 



V 



D ■■ 



4 0.5488(4) 

5 0.3620(3) 

6 0.2292(3) 

7 0.1412(3) 

8 0.0844(4) 

9 0.0505(4) 

10 0.0293(4) 

11 0.0170(3) 

12 0.0097(2) 

13 0.0053(1) 

14 0.00304(6) 

15 0.00179(4) 



0.3166(3) 

0.1420(2) 

0.0593(2) 

0.0233(2) 

0.0087(2) 

0.00315(6) 

0.00111(2) 

0.000380(8) 

0.000128(3) 

5.2(4) x 10~ 5 

1.7(3) x 10~ 5 



0.1888(2) 
0.0591(2) 
0.01648(6) 
0.00424(6) 
0.00101(2) 
0.000226(5) 
5.2(2) x 10~ 5 
1.0(1) x 10~ 5 
2.7(7) x 10- (i 



0.1153(2) 

0.02522(8) 

0.00487(6) 

0.00076(1) 

0.000129(3) 

1.78(7) x 10~ 5 

2.5(4) x 10~ 6 



TABLE V: B k [4, 



k 


D = 2 


D = 3 


D = 4 


D = 5 


4 


-0.01644(5) 


-0.02981(9) 


-0.0370(1) 


-0.0391(1) 


5 


-0.0264(3) 


-0.0316(2) 


-0.0270(2) 


-0.0189(2) 


6 


-0.0285(5) 


-0.0219(4) 


-0.0117(2) 


-0.0059(1) 


7 


-0.0239(5) 


-0.0114(2) 


-0.00403(8) 


-0.00130(3) 


8 


-0.0183(4) 


-0.0056(1) 


-0.00120(5) 




9 


-0.0123(3) 


-0.0025(1) 






10 


-0.0086(4) 








11 


-0.0056(2) 








12 


-0.0038(2) 








13 


-0.0023(3) 









In contrast to the close packed diagrams, we call di- 
agrams Bk[k,i] for which many points are forced to be 
more than a distance of u apart "loose packed" . One se- 
quence of loose packed diagrams is the k point ring of / 
bonds which we denote by R, for which numerical values 
are shown in Table VI. We have found numerically for 
fixed k that as D — > 00 this ring diagram is larger than 
all other studied. We therefore make the following 
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Conjecture: For fixed k, lim^^oo B k /R = 1. In par- 
ticular the sign of B k for large D is (— l) fe_1 . 

For fixed D, we find that for k up to 7 the largest di- 
agrams are of the form of a ring with "insertions" where 
one point of the ring is replaced by a cluster of points con- 
nected by / and / bonds. We thus may define a sequence 
of loose packed diagrams which consist of ring diagrams 
with an inserted diagram. Note that the contribution of 
these sequences to the virial coefficient alternate in sign, 
whereas close packed diagrams always contribute with 
the same sign. The smallest possible insertions have 4 
points, of which there are two examples. The larger of 
these is given in Table VII where the type of insertion is 
indicated in the figure caption. Multiple insertions are 
possible but never dominate for the values of k we have 
studied. 



TABLE VI: B k [k, 1]/-B£ _:i = R/B^ 1 . The underline marks 
the approximate location of the minimum value. 



k 


D = 2 


D = 3 


D = 4 


D = 5 


4 


-0.01639(9) 


-0.0298(1) 


-0.0371(1) 


-0.03925(5) 


5 


0.00860(6) 


0.01623(9) 


0.0214(2) 


0.0230(1) 


6 


-0.00526(8) 


-0.0109(1) 


-0.0150(2) 


-0.01689(9) 


7 


0.00335(6) 


0.0078(2) 


0.0124(2) 


0.0142(3) 


8 


-0.00234(5) 


-0.0064(1) 


-0.0106(2) 


-0.0129(3) 


9 


0.00177(4) 


0.0053(1) 


0.0098(2) 


0.0126(2) 


10 


-0.00125(3) 


-0.00452(9) 


-0.0091(2) 


-0.0126(3) 


11 


0.00095(2) 


0.00392(8) 


0.0089(2) 


0.0128(3) 


12 


-0.00074(1) 


-0.00333(7) 


-0.0083(2) 


-0.0134(3) 


13 


0.00055(1) 


0.00313(8) 


0.0086(2) 


0.0142(3) 


14 


-0.00041(1) 


-0.0027(1) 


-0.0086(2) 


-0.0166(3) 


15 


0.00033(2) 


0.0026(1) 


0.0087(3) 


0.0183(4) 



TABLE VII: B h \k,2\IB\- x = R (-<£>-)/ B*- 1 . Th e underline 
marks the approximate location of the minimum value. 



k 


D = 2 


D = 3 


D = 4 


D = 5 


5 


-0.01017(4) 


-0.01654(7) 


-0.01748(6) 


-0.01551(4) 


6 


0.00756(6) 


0.01191(8) 


0.01257(6) 


0.01067(6) 


7 


-0.00607(8) 


-0.01003(7) 


-0.0106(1) 


-0.00936(9) 


8 


0.0051(1) 


0.0088(2) 


0.0100(1) 


0.00902(8) 


9 


-0.00427(8) 


-0.0083(1) 


-0.0100(1) 


-0.0095(2) 


10 


0.00381(8) 


0.0082(2) 


0.0106(2) 


0.0106(2) 


11 


-0.00309(3) 


-0.0075(1) 


-0.0109(2) 


-0.0116(2) 


12 


0.00259(3) 


0.00728(7) 


0.0115(2) 


0.0136(3) 


13 


-0.00220(3) 


-0.0071(3) 


-0.0125(3) 


-0.0156(3) 


14 


0.00188(4) 


0.0068(4) 


0.0140(5) 


0.0172(4) 


15 


-0.00151(5) 


-0.0060(5) 


-0.0145(7) 


-0.0227(7) 



The data of Tables IV- VII differ in the two cases D = 2 
and D > 3 in the following respect. For D = 2 the close 
packed diagrams Bk[0, 1] and B k [A, 1] arc larger than 
the loose packed diagrams B k [k,i] we have considered, 
whereas for D > 3 the loose packed diagrams rapidly 



dominate the close packed diagrams for large k. Conse- 
quently we will consider the two cases separately. 

For D = 2 an examination of Tables IV- VII indi- 
cates that all diagrams are roughly of the same order 
of magnitude. A measure of the relative size of the close 
packed to the loose packed diagrams is given by the ratio 
\B k [k,l]/B k [0, 1]| which increases from 0.0299 to 0.184 
as k increases from 4 to 15 (0.382 when k = 17). It 
can be argued that if this increase continues then even- 
tually the loose packed diagrams will dominate, but such 
a conclusion requires further computations. The data 
can be interpreted to say there is a qualitative difference 
between D = 2 and D > 3. This interpretation may be 
correct but is not compelling. 

For D > 3 the loose packed diagrams B k [k,i] rapidly 
dominate for large k the close packed diagrams Bk[m, i] 
for fixed m and in particular \B k [k, 1]| first becomes 
larger than B k [0, 1] for k = 9 in D = 3, k = 7 in D = 4, 
and k — 6 in D — 5. Further studies indicate that B k [k, i] 
even dominates B k [k — 1, i) for large k. Thus in order to 
determine the radius of convergence we will restrict our 
attention to the diagrams belonging to B k [k, i}. 

We see in Table VI that the Ree-Hoover ring diagram 
B k [k, lJ/Bg -1 increase in magnitude for D > 4 when k 
is sufficiently large. In D = 3 the data of Table VI can 
be interpreted as cither having a minimum near k = 1 or 
possibly approaching a constant as k — > oo. 

The largest diagram shown is the four point insertion 
B k [k, 2] of Table VII. It seems clear that for D > 4 this 
diagram will always be bigger than the Rcc- Hoover ring 
but the ratio seems to grow only as a power of k. We 
believe the same behavior happens for D = 3 but it is 
more difficult to see with the precision given in the tables. 

From these and other numerical studies we conjecture 
that the exponential rate of growth of the ring and ring 
diagrams with insertions is the same. 

To estimate the rate of growth we concentrate on the 
Ree-Hoover ring of Table VI. We note that the absolute 
value of this diagram must be strictly less than the ab- 
solute value of the Mayer ring diagram which is obtained 
by replacing all the / bonds by unity. The Mayer ring 
may be expressed as a single integral [41], and from this 
we obtain 



B k [k,l] 



< 



{k- 1)(2tt)^ 



dxx 



D-l 



Jp/2{x) 

x D / 2 



(3) 

where = 2tt d / 2 /T(D/2) and J D/2 (k) is the Bessel 

function of the first kind. The large k behavior of this 
integral is obtained by steepest descents and thus we have 

as k — > oo 



B k [k,l] 



B. 



fe-i 



< 



(fc- l)(l + £>/2) D / 2 
fe 1 + D / 2 r(l + D/2) 



,fc-2 



(4) 



From Table VI we see that for all D the ratios B k+ i [k+ 
1, l]/(B 2 B k [k, 1]) are all substantially below this bound 
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of 2. If the data in Table VI are extrapolated to say for 
D = 3 that Bk[k, 1] goes to a constant (or a power of k) 
as k — ► oo then we see that the radius of convergence of 
the sum of Ree-Hoover ring diagrams is r\ r h = 0.25. At 
worst the ratios in Table VI for D = 3 are bounded below 
by 0.91 as k — > oo which leads to a radius of convergence 
of f] T h = 0.27. This radius of convergence is substantially 
less than the freezing density of r/f = 0.49. Similarly we 
estimate from Table VI that in D = 4 the radius of con- 
vergence of Ree-Hoover rings is r\ r h ~ 0.12 and in D = 5 
is i] r h ~ 0.052 which are to be compared with the freezing 
densities rjf = 0.31 in D = 4 and i]f = 0.19 in D = 5 ob- 
tained from [42-44]. The alternation of sign of the ring 
means that the leading singularity of the sum of these 
diagrams is on the negative r\ axis. 

If all loose packed diagrams Bk[k,i] had the sign 
(— l) fc_1 of the Ree-Hoover ring diagram then the esti- 
mates of r] r h found above would be an upper bound on 
the radius of convergence of the virial expansion. How- 
ever for every order k both signs occur in the class Bk [k, i] 
and in particular the diagram Bk[k, 2] is larger than and 
has the opposite sign from Bk[k,l}. Therefore cancella- 
tions of diagrams within the class Bk[k,i] can occur and 
as evidence of such cancellation we note that the virial 
coefficient B%jB\ of Table I for D = 3 is smaller than 
both B 8 [8, 1] and Bg[8,2], and has the opposite sign to 
the ring diagram B$[8, 1]. 

However, the diagrams Bk[k,2] also alternate in sign, 
and by themselves lead to a singularity on the negative 
real axis. If this singularity is to be avoided extensive 
cancellation must take place beyond what can be seen 
in the virial coefficients up through order 8 as given in 
Table I. From this point of view we see that a detailed 
study of diagrams for orders substantially greater than 
8 is needed to substantiate any claim concerning the ra- 
dius of convergence of the virial series for hard spheres. 
None of the approximate equations of state for hard 
spheres [3, 14, 18, 21-30] includes diagrams of such high 
order. We therefore conclude that there is no existing 
evidence to support the claim that the virial expansion 
converges beyond the freezing density ijf of hard spheres. 

This work was supported in part by the National 
Science Foundation under DMR-0073058. We thank 
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* Electronic address: Nathcm.Clisby@stonybrook.edu 
^ Electronic address: mccoy@insti.physics.sunysb.edu 
[1] J. L. Lebowitz and O. Penrose, J. Math. Phys. 5, 841 
(1964). 

[2] B. J. Alder and T. E. Wainwright, J. Chem. Phys. 27, 
1208 (1957). 

[3] W. G. Hoover and F. H. Ree, J. Chem. Phys. 49, 3609 
(1968). 



[4] F. H. Ree and W. G. Hoover, J. Chem. Phys. 40, 2048 
(1964). 

[5] J. G. Loeser, Z. Zhen, S. Kais, and D. R. Herschbach, J. 

Phys. A. 26, 4805 (1993). 
[6] M. Luban and A. Baram, J. Chem. Phys. 76, 3233 

(1982). 

[7] M. Bishop, A. Masters, and J. H. R. Clarke, J. Chem. 

Phys. 110, 11449 (1999). 
[8] L. Tonks, Phys. Rev. 50, 955 (1936). 
[9] L. Boltzmann, Verslag. Gewonee Vergadering Afd. Natu- 
urk. Nederlandse Akad. Wtensch. 7, 484 (1899). 
[10] J. S. Rowlinson, Mol. Phys. 7, 593 (1964). 
[11] P. C. Hemmer, J. Chem. Phys. 42, 1116 (1964). 
[12] J. J. van Laar, Amsterdam Akad. Versl. 7, 350 (1899). 
[13] B. R. A. Nijboer and L. van Hove, Phys. Rev. 85, 777 
(1952). 

[14] F. H. Ree and W. G. Hoover, J. Chem. Phys. 40, 939 
(1964). 

[15] K. W. Kratky, J. Stat. Phys. 27, 533 (1982). 
[16] K. W. Kratky, Physica A 87, 584 (1977). 
[17] K. W. Kratky, J. Stat. Phys. 29, 129 (1982). 
[18] E. J. Janse van Rensburg, J. Phys. A 26, 4805 (1993). 
[19] F. H. Ree and W. G. Hoover, J. Chem. Phys. 46, 4181 
(1967). 

[20] A. Y. Vlasov, X.-M. You, and A. J. Masters, Mol. Phys. 

100, 3313 (2002). 
[21] E. Thiele, J. Chem. Phys. 39, 474 (1963). 
[22] M. S. Wertheim, Phys. Rev. Lett. 10, 321 (1963). 
[23] E. A. Guggenheim, Mol. Phys. 9, 199 (1965). 
[24] N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 

635 (1969). 

[25] R. Hoste and J. D. Dael, J. Chem Soc. Faraday Trans. 2 
80, 477 (1984). 

[26] J. I. Goldman and J. A. White, J. Chem. Phys. 89, 6403 
(1988). 

[27] E. J. L. Fevre, Nature Phys. 235, 20 (1972). 

[28] D. Ma and G. Ahmadi, J. Chem. Phys. 84, 3449 (1986). 

[29] S. Jasty, M. Al-Naghy, and M. de Llano, Phys. Rev. A 

35, 1376 (1987). 
[30] Y. Song, R. M. Stratt, and A. E. Mason, J. Chem. Phys. 

88, 1126 (1988). 
[31] J. D. Bernal, Proc. Roy. Soc. Lond. A 280, 299 (1964). 
[32] J. G. Kirkwood and E. Monroe, J. Chem. Phys. 8, 845 

(1940). 

[33] T. V. Ramakrishnan and M. Yussouff, Phys. Rev. B 19, 

2775 (1979). 
[34] A. D. J. Haymet and D. W. Oxtoby, J. 

1769 (1986). 
[35] H. N. V. Temperley, Proc. Phys. Soc. 

536 (1957). 

[36] W. G. Hoover and A. G. D. Rocco, J. Chem. Phys. 36, 
3141 (1962). 

[37] R. J. Baxter, J. Math. Phys. 11, 3116 (1970). 
[38] R. J. Baxter, J. Phys. A 13, L61 (1980). 
[39] G. S. Joyce, Phil. Trans. R. Soc. Lond. A 325, 643 (1988). 
[40] F. H. Ree and W. G. Hoover, J. Chem. Phys. 41, 1635 
(1964). 

[41] S. Katsura and Y. Abe, J. Chem. Phys. 39, 2068 (1963). 
[42] J. P. J. Michels and N. J. Trappeniers, Phys. Lett. 104, 

[43; 

[44 



Chem. Phys. 84, 
(London) B 70, 



S. Katsura and Y. Abe, 
J. P. J. Michels and N. 
425 (1984). 

M. Luban and J. P. J. 
(1990). 

R. Finken, M. Schmidt, 
016108 (2001). 



Michels, Phys. Rev. A 41, 6796 
and H. Lowen, Phys. Rev. E 65, 



